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Abstract. Given the extreme accuracy reached by future global space astrometry, one 
needs a global relativistic modeling of observations. A relativistic definition of astrometric 
observables is then essential to find uniquely coordinates, parallax and proper motion of a 
stellar object. The standard procedure is to explicitly solve the null geodesic equations to 
describe the trajectory of a photon emitted by a celestial object and received by a moving 
observing satellite. However, this task can be avoided if one builds an astrometric set up 
using the recent formalism of time transfer functions. 
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1. Introduction 

Future space astrometry missions, such as 



Gaia (Bienayme & Turon 2002|l, will soon pro 



vide large astrometric catalogues with typi- 
cal accuracy of several microarcseconds (yuas). 
It is nowadays well-known that //as astrom- 
etry needs a precise relativistic modeling to 
define and calculate accurate astrometric pa- 
rameters. In particular, when a satellite ob- 
serves a celestial object, one needs to know 
its attitude and onboard relativistic proper time 
as well as the deflection of the incoming 
light ray detected on the CCD. In the context 
of the Gaia mission, several approaches are 
available, notably GREM (|Khoner|[2004] l and 
RAMOD flde Felice et al.|2006) . Both model- 
ings need to determine the full light trajectory 
from the emitting celestial object to the observ- 
ing satellite by solving the null geodesic equa- 
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tions. However, it has been recently demon- 
strated that this task is not at all mandatory 
and can be replaced by another approach, ini- 
tially based on the Synge World Function 
( |Le Poncin-Lafitte et al.||2004|> and then on 



the time transfer functions (Teyssandier & 
Le Poncin-Lafitte 2008). The purpose of this 



article is to build straightforwardly an astro- 
metric set up using the formalism of the time 
transfer functions. 



This paper is organized as follows. 
Section [2] gives the notations used in this arti- 
cle. In sectionRlwe set up the astrometric prob- 
lem. Section |4jrecalls the relation between the 
covariant components of the vector tangent to 
a light signal and its flight time between two 
events. This relation is then explicitly given in 
section [5] within the post-Newtonian approxi- 
mation. Finally, section [6] delivers some con- 
cluding remarks. 
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2. Notations and conventions 

Throughout this work, c is the speed of light in 
vacuum and G is the Newtonian gravitational 
constant. The Lorentzian metric of space-time 
V4 is denoted by g. We adopt the signature 
(H ). We suppose that space-time is cov- 
ered by some global coordinate system x a = 
(x°,x), with x° = ct and x = (x 1 ), centered 
on the Solar System barycenter. Greek indices 
run from to 3 and Latin indices from 1 to 

3. Moreover, we assume that the curves of 
equation x' = const are time-like, at least in 
the neighborhood of the chosen observer. This 
condition means that goo > in the vicin- 
ity of the observer. Any ordered triple is de- 
noted by a bold letter. In order to distinguish 
the triples built with the space-like contravari- 
ant components of a vector from the ones built 
with covariant components, we systematically 
use the notation a = (a 1 , a 2 , a 3 ) = (a 1 ) and 

b = {b l ,bM = Q> i ). 



3. Set up of the astrometric problem 

Let us consider a time-like observer. All along 
its worldline, we choose a tetrad of four vec- 
tors E'fy where index (a) is only the tetrad in- 
dex, which runs from to 3 to enumerate the 4- 
vectors, while fi is a normal tensor index which 
can be lowered and raised by the use of the 
metric tensor : 



orthogonal to is given by Teyssandier & 



|Le Poncin-Lafitte 



[((2006} 



k!' 



F v k 



£ (0)- 



(3) 



Then, we calculate the three director cosines 
cos <p( a ) formed by each space-like vector EY-. 

(here index (a) runs from 1 to 3) with W 



cos ( 



k^E, 



(a) ■ 



(4) 



Taking Eqs. (|3]-[4]) into account, the final ex- 
pression of cos 0( fl ) can be written as follows 



cos ( 



Ei , "H ki Ei / , 

(a) ' (a) 



E° 

£ (0) 



(5) 



where /3' = v'/c, v' is the coordinate velocity of 
the observer and k t = kj/ko will be called in the 
following the deflection functions. 

Eq. Q shows that we only need the deflec- 
tion functions at the reception event to deter- 
mine the director cosines of an astrometric ob- 
servation and that it's not necessary to trace the 
light ray from the emission event to the recep- 
tion event. In the next session, we recall the re- 
lation between time delay and light deflection. 



E ( a ) M - g^v E v (a) , E^ a) - g^E^y. ( 1 ) 

The tetrad vectors are required to satisfy the 
following property 

which implies that the vectors are orthogonal 
to each other, vector £^ 0) being unit and time- 
like, and E*L being unit and space-like. In the 
following, we do not specify more additional 
information about this tetrad. 

Let us now study the reception of a light ray 
from a celestial object by a satellite. Its direc- 
tion is given by the vector kf = (k°, k'), tangent 
to the incoming light ray and its unit space- 
like direction k^ relative to the hyper-surface 



4. From coordinate time delay to light 
deflection 

Let us consider that space-time is globally reg- 
ular with the topology MxM 3 and that it is with- 
out horizon, which is admissible in the context 
of practical space astrometry within the Solar 
System. We suppose the existence of a unique 
light ray connecting xa = (cf^,^), the event 
of emission of a photon from a celestial object, 
and xb = (cfg,^), the event of reception of 
the light by the observing satellite. The differ- 
ence ts - tA is the (coordinate) travel time of 
the photon between the two points. This quan- 
tity may be considered either as a function of 
the instant of reception tg and of Xa, Xb, or as 
a function of the instant of emission tA and of 
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x A and x B . So, it is possible to define two time 
transfer functions T r and T e by putting] 

t B -t A = T r (x A , t B , x B ) = T e (t A ,x A ,x B ), (6) 

where T, and T e are called the reception 
and the emission time transfer function, re- 
spectively. Explicit expressions of these func- 
tions are different except in a stationary space- 
time in which the coordinate system is chosen 
so that the metric does not depend on x°. Here, 
we are obviously most interested in T r . Indeed, 
when the satellite is doing an observation, we 
know its state's vector and its onboard time. 

Moreover, a theorem relating the de- 
flec tion functions and T, has been found 
by (Le Poncin-Lafitte et al. |2004| l, where fc is 
obtained from the partial derivatives of T r as 
follow 



(k.) =-, 



dx' 



1 - 



~dt~B 



(7) 



The question is then how to determine this 
function in order to calculate the astrometric 
director cosines. 

5. Post-Newtonian time transfer and 
deflections functions 

It was shown that T r is the solution 
of the following partial differential equa- 
tion ( Teyssandier & Le Poncin-Lafitte|2008 ) 



g™(x B - cT r ,x A ) + 2cg m (x B - cT r ,x A ) 

A, dx' A d^ 



9T r 
dX 



A 

(8) 



Finding a solution to this equation is 
as challenging than a determination of the 



null geodesic equations (Kopeikin & Schafer 
1999 ). However, in the weak field approxima- 
tion, this task is easier. Let us write the metric 
tensor as follow 



(9) 



1 Note that the order of the arguments of T r 
in E q. ((o} is different from th e one given in 



Ref. ( |Le Poncin-Lafitte et al.|2004| . 



with 77 a Minkowskian background and h a per- 
turbation admitting the following general post- 
Minkowskian expansion 



(10) 



It was possible to show that T r can be writ- 
ten as a series of ascending powers of the 
Newtonian gravitational constant G 

R °° 

T r (t B ,x A ,x B ) = - + Y \G n T} n \t B ,x A ,x B ),(U) 

where it has been demonstrated that each per- 
turbation term 7/ can be obtained as a line 
integral along a null straight line. This result is 
particularly interesting and can be interpreted 
as a Fermat principle (Perlick 1990) in the «th- 
post-Minkowskian approximation. 

In particular, if we work in the con- 
text of the post-Newtonian approximation, 
and taking the metric tensor recommended 
by the International Astronomical Union in 
2000 ( |Soffel et al.|2003| ), we can write the met- 
ric tensor g„ v for a light signal as 



*oo = l-- i W(t,x)+0(c- 4 ) , 
g 0i = -^W i (t,x)+O(c- 5 ) , 



(12) 



where y is a PPN parameter (Will 1993 1, 
W(t,x) and W = {W(t,x)} being the grav- 
itational potentials at the coordinate (f, x). 
After substitution of the metric into Eq. (|8]), a 
straightforward calculation leads to 



TritsiXA^E) = — 
C 



r r l 

"?Jo 



(y + W- -W-N 



(13) 

dA + 0(c- 5 ), 



where the integral is taken along the null 
straight line 



Z-(X) = (ct B — AR,x B — ARN) ., 



(14) 
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with/? = \x B -x A \ and A? = = (x B -x A )/R 
being the Minkowskian direction between x A 
and Xb- 



Considering Eqs. ( 13 - 14 1, it is possible to 
calculate the partial derivatives of T r . One gets 



d, B T r 



(y+l)R 



f 

Jo 



8W 



at 



dA 



(15) 



4R 



and 



AT 1 



(y+1) 



4 f> 

c 4 Jo 



lVAf ! ' + (l -A)R 



dW 

~dx l 



(16) 



f (r+i) ,<3W 

W + V/ - J ARN' — 
4 <9r 



+(1 -/I)/? A 7 



3W 



^+(9(c" 5 ). 



Substituting Eqs. ( 15 - 16 1 into Eq. Q, we ob- 
tain finally the following relation for the de- 
flection functions 



-N' 



(y + i) 



(17) 



I' 

Jo 



quite new and needs further development but, 
in principle, it avoids the difficult task of solv- 
ing the null geodesic equations. 
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WN' + (1 - A)R 



dW 



dA 



i 

Jo 



(y+1) dW 

W-- -(l-A)RN'^- 

4 ot 



dA + 0(c- A ). 



6. Conclusions 

We presented in this article a relativistic pro- 
cedure for global astrometry based on the time 
transfer functions formalism. We were able to 
obtain the light deflection of a light ray con- 
necting a stellar source and an observing satel- 
lite solving the boundary value problem. The 
deflection functions in the Post-Newtonian ap- 
proximation were then obtained as a line in- 
tegral of the gravitational potentials along the 
Minkowskian line of sight. This approach is 



